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Abstract. The two most common mass estimators that 
use velocity dispersions and positions of objects in a stel- 
lar system, the virial mass estimator (VME) and the pro- 
jected mass estimator (PME), are revisited and tested us- 
ing A^-body experiments. We consider here only spherical, 
isolated, and isotropic velocity dispersion systems. 

We have found that the PME can overestimate masses 
by « 20%, for realistic cluster mass profiles, if applied only 
to regions around the total system's effective radius. The 
VME can yield a correct mass at different radii provided 
that the potential energy term is correctly taken into ac- 
count and the system is completely sampled; otherwise, it 
may lead to similar errors as the PME. 

A surface pressure {3PV) term recently alluded to be 
usually neglected when using the VME and therefore re- 
quired as a correction term is here found not necessarily 
required, although it can be used to yield a reasonable 
correction term. The preferred method here, however, is 
the virial theorem due to its simplicity and better agree- 
ment with A-body experiments. The possible reasons for 
the mass discrepancies found when using the PME and 
the VME in some A-body simulations are also briefly dis- 
cussed. 
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1. Introduction 

Clusters of galaxies, probably the largest known virialized 
systems, provide an important tool to study the large-scale 
structure of the universe (e.g. Peebles 1980); for example, 
they provide an estimate for the dynamical mass density 
parameter on scales of ~ 1 Mpc. Clusters can also distort 
the images of faint background galaxies by weak lensing 
providing a mean to study the cluster mass distribution 
and thus placing constraints on large-scale structure for- 
mation models (e.g. Miralda-Escude 1995). 
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Several methods are used to determine the mass of 
clusters of galaxies: e.g. velocity dispersions (Bahcall 
1977); X-rays of the intracluster gas (Jones & Forman 
1984); and gravitational lensing (Tyson et al. 1990). In 
general, although the three previous methods yield similar 
results for the total mass of rich clusters, ^ 1O"'^^~^^M0, 
they provide discordant values for the mass at different 
radii. This provides a motivation to have well character- 
ized and tested methods for total mass and mass profile 
determination, owing to their importance e.g. in cosmol- 
ogy (White et al. 1993). We will deal here only with the 
first method mentioned. 

Analysis of the structural parameters in clusters of 
galaxies show that they are not really spherically symmet- 
ric structures, but their shape appears to be somewhat 
triaxial (e.g. Plionis et al. 1991). Moreover, since some 
clusters still appear to be on the process of collapse, e.g. 
judging from the presence of substructures (West 1994), 
their velocity distribution is probably not isotropic. This 
may be particularly true at their outer parts where galax- 
ies or material can be still infalling with predominantly 
radial orbits (White 1992). 

Nevertheless, observational limitations and simplicity 
makes us to assume virial equilibrium and sphericity for 
clusters when estimating their mass or mass profile. But 
even in such ideal scenario, methods that use observed ra- 
dial velocities and projected positions are still a matter 
of some discussion. On the one hand, for example, Heisler 
et al. (1985, hereafter HTB) concluded that the projected 
mass estimator (PME) is probably the best option for es- 
timating the mass of clusters of galaxies from galaxy mo- 
tions. On the other hand, Perea et al. (1990) reached the 
conclusion that the best estimator is the virial mass esti- 
mator (VME) and that, for example, anisotropy, a mass- 
spectrum, or the presence of substructure can lead to an 
overestimation of the mass if they are not properly taken 
into account. 

Furthermore, it seems that the above mass estimators 
do not even provide adequate masses when applied to A- 
body simulations; this is more critical, since in these sys- 
tems we have complete information on the positions and 
velocities of particles. For example, Thomas & Couchman 
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(1992, hereafter TC) using A^-body simulations of the for- 
mation of clusters of galaxies report that mass estimates 
based on the virial theorem underestimate the total mass 
by a factor w 4, and that the PME yields a factor w 2 too 
small. TC conclude that the PME is probably the best 
option to determine the mass using velocity dispersions. 

More recently Carlberg et al. (1997a, hereafter CYE) 
and Carlberg et al. (1996) have raised concerns regarding 
the consistency of the usual form of the VME when the to- 
tal system is not entirely sampled. In their investigation of 
the average mass of galaxy clusters they establish that the 
VME overestimates their total mass by « 20%, attribut- 
ing this discrepancy to the neglect of the surface pressure 
(3PV) term in the continuous form of the virial theorem. 
A correction to virial mass estimates based on this 3PV 
term is starting to be applied to clusters of galaxies by 
the community (e.g. Girardi et al. 1998). Note that this 
finding of CYE is contrary to the trend observed by TC 
for the virial mass estimate. 

From our investigation of the VME and the PME pre- 
sented below, we will show that although the 3PV term 
can account for the mass overestimation when the VME 
is applied to a subsample of an equilibrium gravitational 
system, the correct application of the VME yields also a 
correct mass at different radii. The physical reason for the 
overestimation of mass by the VME in the previous situ- 
ation lies in an incomplete consideration of the potential 
energy of the subsystem. We will show that when the po- 
tential energy is fully accounted for no discrepancy exists 
at any radius between the VME and the true mass for an 
A'^-body system. In an astronomical situation, however, 
the previous result relies, of course, on having knowledge 
of the total extent of the system which is a somewhat am- 
biguous matter; this relates to the problem raised by CYE. 
In practical terms, this also affects the use of the 3PV 
term as a correction factor since one needs a fair knowl- 
edge of the system's equilibrium extent, in addition to the 
velocity dispersion profile of a cluster. On other hand, the 
VME has the virtue of not depending on the form of the 
orbital distribution of galaxies in a near-equilibrium clus- 
ter, hence being straightforward to apply. 

Haller & Melia (1996) have considered also the prob- 
lem of calculating the mass profile by use of the PME. One 
of their conclusions is that a boundary term, arising from 
the finite sampling of a gravitational system, may make 
appreciable contributions to the mass estimate at the in- 
ner regions of stellar systems. Their conclusion is similar 
to the one we reach here in an independent manner, but 
this boundary term is studied here to a larger extent using 
different 'cluster' models and the physical reason for the 
overestimate of the mass is elucidated. 

In this work we revisit and test the PME and VME 
for an isolated spherical system of identical particles with 
an isotropic velocity distribution. We will quantify in par- 
ticular the effect of a boundary term usually neglected in 
the PME and provide bounds to the errors that one may 



have when using the VME on a subset of an equilibrium 
system with different profiles. In a future paper we will 
apply our results to obtain mass profiles of nearby clus- 
ters and compare them with e.g. the profiles derived from 
X-rays. 

The plan of the paper is as follows. In §//, the stan- 
dard PME formalism is revisited and tested using different 
spherical models considered appropriate for the descrip- 
tion of clusters of galaxies; e.g. those having a 'cusp' and a 
'core'. A correction term is explicitly obtained and tested 
with A'^-body experiments. In §///, the virial theorem, 
both in its discrete and continuous form, and the effect of 
the 3PV term on the mass determination are investigated. 
Finally, in ^IV, our main conclusions are presented. 

2. The Projected Mass Estimator 

This estimator uses the average of the quantity v^R over 
the whole tracer population (e.g. galaxies in a cluster) 
relating it to the total mass by (HTB, Perea et al. 1990) 
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where Vzi is the observed line-of-sight velocity of a galaxy 
relative to the cluster mean, and Ri is its projected radius 
from the centre of the distribution. The factor fp is re- 
lated to the orbital distribution of galaxies and therefore 
a function of the anisotropy parameter (3 (see below), be- 
ing equal respectively to 64/7r and 32/7r for the cases of 
radial {(3 = 1) and isotropic {(3 = 0) orbits. 

Equation (1) is obtained as follows. On the one hand, 
for a spherically symmetric system under steady-state con- 
ditions and no rotation the Jeans equation holds (Binney 
& Tremaine 1987) 
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where p is the mass density, ip the gravitational potential, 
(Tr the radial velocity dispersion, and (3 — I — ^ll'^l the 
anisotropy function; at is the tangential velocity disper- 
sion. Multiplying Jeans equation by r^, and integrating 
by parts up to a radius r we have: 
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On the other hand, the average {v^R) in Eq. (1) up to 
a particular radius r is estimated as (e.g. HTB): 



{vlR), 
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' 9 -\- a'^ sin^ 9 is the velocity dispersion 
along the line-of-sight, R = r sin 9, and /(r, v) is the phase 
space density, defined such that /(r, v)d'^rd'^v is the mass 
contained in the phase space volume element d'^rd'^v. 
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Evaluating the integral in Eq. (4), by using Eq. (3) and 
assuming isotropy (/3 = 0), we have: 



{<R)r = 



4M(r) 



r'^p(r)a^(r) + I r^p(r)— — dr 
dr 



(5) 



Since d(p = GM{r)dr/r'^ and Airr'^p = dM/dr we have, 
after rearranging terms: 

A(r) EE ^{vlR)r - M{r) = ^^ r'p{r)cTl{r) . (6) 

This function A(r) gives us directly the effect of neglecting 
the boundary term on the estimation of the mass via the 
PME when partial sampling is done. Now, since one usu- 
ally fits a particular projected model to a set of data, the 
quantities to evaluate A(r) follow from the correspond- 
ing deprojected model. One may use e.g. Eq. (9) to eval- 
uate the radial velocity dispersion in case this does not 
exists explicitly beforehand for the adopted model. Inte- 
grating over the whole system we recover the standard 
formula of the projected mass for an isotropic stellar sys- 
tem: Mp = 32(t;2i?)/(7rG). 

Therefore, in using equation (1) for an isotropic system 
one must remember that this form of Mp holds only when 
the total extent of the stellar system has been sampled; 
otherwise, the correction term A(r) has to be applied. As 
will be shown below, in astronomical situations the use 
of Mp can yield a significant overestimate of the system's 
mass especially if we sample only around its effective ra- 
dius. Furthermore, since the boundary of an astronomi- 
cal system is somewhat ambiguous, particularly if large 
amounts of dark matter exist at the outer optical edges, it 
is important to quantify the error one might commit when 
using the PME. We proceed to estimate such error in the 
next subsection. 

2.1. Boundary Term for Spherical Systems 

In this subsection we estimate theoretically the boundary 
term A(r) for four different spherical models: (1) de Vau- 
couleurs' or R^/'^ (de Vaucouleurs 1948), (2) Hernquist's 
(1990), (3) Dehnen's (7 == 0) (1993), and (4) a King's mod- 
ified profile. We consider these models to bracket most of 
the possible profiles that can be fitted to a cluster of galax- 
ies, which will be here the astronomical object of interest. 
Evidently, these profiles can be applied to other spherical 
systems, as well as the mass estimators. 

Navarro, Frenk & White (1996), hereafter NEW, have 
proposed a profile model appropriate for clusters of galax- 
ies of the form p{r) ex r^^ {r + a)^"^ , with a being a scale 
radius. Carlberg et al. (1997b) have shown that the av- 
erage mass density profile of clusters is well described by 
the NEW profile though they mention that the Hernquist's 
model is statistically acceptable. However discordant re- 
sults exist in the literature on the subject of cluster profiles 
(e.g. Adami et al. 1998). We have not used the NEW pro- 
file here due to its intrinsic difficulty in defining formally 



a total mass, hence a half-mass or effective radius, unlike 
the other profiles. Indeed, in the NEW profile p -^ r~^ 
as r ^ 00 and the total mass diverges. We do not ex- 
pect much difference in the boundary term A(r) when 
the NEW and Hernquist profiles are fitted to real clus- 
ters. Thus, for convenience, we deal here only with the 
Hernquist's model. 

Here the units adopted for the numerical calculations 
are such that the total mass is unity, M = 1, the total 
system's effective radius i?o = 1, and the gravitational 
constant G = 1. Following, for completeness, we write 
some pertinent quantities to compute A(r) for the differ- 
ent models to be considered. 

2.1.1. De Vaucouleurs Model 

The surface brightness of the de Vaucouleurs' model is: 

/(i?)-/eexp{-M(i?/i?e)^/'-l]}, (7) 

where Re is the effective radius, a scale-radius where half 
of the total light is emitted, and b — 7.66925. We will 
consider here the effective radius to be the locus where half 
of the projected mass resides. The total mass is, assuming 
a mass-to-light ratio of unity, M ~ 8\tt cxp(b)R^Ic/b^ ~ 
22.7i?2/^. 

The mass density p{r) can be obtained by deprojecting 
the surface brightness using the Abell integral equation 
(Binney & Tremaine 1987). The latter has been obtained 
in implicit form by Poveda et al. (1960) as: 
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where t = i?^/^, and 

The product pa^ can be obtained from the Jeans equa- 
tion by imposing the condition that pa^. ^ as r ^ 00. 
Since p and dip/dr are non-negative, and a^. always posi- 
tive, we have: 



dpo-^ 
dr 



= -P-^ -^ P{r)<^r{r) = / 



^(^dr. (9) 



The mass M{r) follows from integrating the mass density, 
and the boundary term can then be computed numerically. 

2.1.2. Hernquist Model 

This model is similar to an i?^' * model, but has the virtue 
of being expressed in terms of simple functions. Some of 
the quantities of interest here are: 
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where a is a scale-radius, and M the total mass. The ef- 
fective radius here is Re = 1.8153a . 

Using Hernquist's expressions and integrating several 
times by parts equation (9), we find that: 
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2.1.3. Dehnen 7 == Model 



Dehnen (1993) has provided a family of potential-density 
models for spherical systems. We use his 7 = model 
since it provides a core resembling King models and allows 
simple analytic calculations (see also Tremaine et al. 1994, 
model 77 = 3). The density and mass are given respectively 
by: 
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(12) 
2/30.1 



The velocity dispersion is ct^ — GM{a + 6r) (r + a 
The effective radius in this model is _Rc = 2.9036a . 

2.1.4. King Modified Model 

The simplicity of this empirical model, besides allowing 
larger cores, makes it well suited for fitting observational 
data. The projected and physical density, and projected 
mass are (e.g. Perea et al. 1990, Adami et al. 1998), re- 
spectively: 
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where a, and Rq, the structural length of the model, are 
fitting parameters. The effective radius, in terms of Rq, 
is i?2 = Rl [2i/("-i) _ 1]. The total mass M may be es- 
timated e.g. by the virial theorem or considered another 
fitting parameter. The following relations among the dif- 
ferent parameters hold: 



M = 
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a — 1 



■ 7r/ni?, 
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y/TTT{a) Ro 



r being the Gamma function. Hereafter we will consider 
only the value of a = 1.3, which is close to fittings to 
the Coma cluster profile (e.g. Perea et al. 1990) and to 
the spectroscopic value of a /?-model of X-ray emission in 
clusters (Sarazin 1988). We refer hereafter to this model 
as King. 

^ A misprint in Dehnen's Eq. A3 is corrected here. 
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Fig. 1. Boundary term in the PME as a function of 
radius for de Vaucouleurs', Hernquist's, Dehnen's, and 
King's models. The maximum values occur respectively at 
^ (0.17,0.18,0.19,0.21) and inside R^ . These values give an 
indication of the error that one may expect to make when es- 
timating masses via the PME at different radii of e.g. a cluster 
of galaxies. 



2.1.5. Theoretical Value of the Boundary Term 

In Fig. 1 we show the function A(r) for the above consid- 
ered models. Of the models considered here the Dehnen's 
model lends itself to an easier computation of the max- 
imum value of the boundary term. Indeed, the following 
expression for the boundary term is readily obtained 



A(r) = 



GMar{a + 6r) 
5 (r + a)3 



(15) 



The maximum of A(r) occurs at Tmax = (5 -I- v31)a/6, 
or ss 0.61 for our chosen units, and has a value of 0.194 . 
For the de Vaucouleurs' model the maximum value of this 
boundary term is 0.166, for the Hernquist's one 0.176, 
and for the King's 0.21 . The radius at which these three 
maxima occur are at 1.07, 0.66, and 0.43, respectively; i.e. 
all occur at r « i?e • 

A trend for the maximum of A(r) to increase as the 
system gets less concentrated is observed. The latter is 
explained if we consider the extreme case of an isothermal 
sphere, where p oc CT^r"^^ M ex cr^r, and a^ is constant, 
from where we obtain that A(r) ex r. Therefore Amax 
grows as the 'isothermal' region of the model grows. Since 
the boundary term is positive and linear in M we can be 
overestimating the total mass by about 20% by using the 
quantity 32{v'^R)/{ttG) as in Eq. (1, /3 = 0), if the sam- 
pling is done only up to regions near the effective radius 
of the gravitational system. The precise value of the error 
depends on the density profile of the system. 
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2.2. Numerical Test of the PME 

Here we evaluate Mp — 32{v'^R)/{ttG) as a mass estimator 
by using ./V-body models. We only present here for conve- 
nience the results for the Dehnen's model, but the same 
quantities and analogous results consistent with the theo- 
retical expectations were found for the other three models 
of §2.1. 

A Monte Carlo realization of A^ = 10^ particles, each of 
mass rrii — 1/N , was constructed. The numerical 'cluster' 
was first evolved using a TreeCode (Barnes & Hut 1986) 
in order to test its stability before the PME is applied 
to the initial configuration. The total time of evolution 
lasted 160 time-steps, each of At ~ 0.031 units, with soft- 
ening parameter e = 0.025 and tolerance parameter 9 = 1. 
Quadrupole terms were included. Energy changed by less 
than 0.1% throughout the evolution. 

The ratio of kinetic to potential energy, as provided 
by the code, was initially T/U = —0.52 reaching T/U = 
—0.50 at the end of the computation. The initial discrep- 
ancy of the numerical model from the ideal virial ratio is 
attributed to the initial positions and velocities not being 
in equilibrium with the code's potential which depends on 
the softening introduced: i.e. iy9 ex {r + e)~^. After this test 
was conducted, we considered our initial Af-body system 
to be in stable equilibrium so we could apply to it the 
mass estimators considered in this work. 

To compute 32{vIR)/{ttG) we proceeded as one would 
do observationally to compute the mass at different 
radii of clusters. We divided the A^-body cluster into 40 
projected concentric radii, spaced logarithmically, from 
(0.1 -^ 100)i?o • Inside each radius the summation of 
vl.Ri, with Ri = {xf + yfY/"^, was calculated for all per- 
tinent particles and then divided by their number N{Ri). 
In this way we estimated 32{v'^R)/{ttG) or, equivalently, 
the observational projected mass Mp. 

In Fig. 2 we plot Mp as a function of projected ra- 
dius for our Dehnen numerical cluster. Also shown are: the 
virial mass (My), as discussed below, the true projected 
mass M{R), and the real 3-D mass within radius r, M{r). 
In the upper panel of Fig. 2, the difference Mp — M{r) is 
shown as a long-dashed line, the difference Mp — M{R) as 
a dotted line, and the theoretical value of A(r) as in equa- 
tion (15) is shown as a solid line. The difference between 
Mp and M{r) can be said to be in perfect agreement, both 
in magnitude and in position of the maximum, with the 
theoretical expectations of §2.1. Poisson noise can explain 
the obtained differences as we have verified, although it 
is not shown in the figure to avoid overcrowding. The dif- 
ference Mp — M{R), although similar in shape to A(r), 
is about 50% smaller in magnitude. A similar behaviour 
of these discrepancies was observed for the other models 
considered here. 
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Fig. 2. Results of the application of the PME (Mp) and VME 
(Mv) as a function of radius to Dehnen's (7 = 0) numerical 
model of a cluster. Also, the theoretical projected mass profile 
M{R) and physical mass M{r) are shown for comparison. The 
upper panel shows the differences Mp — M{r) {long- dashed), 
Mp — M{R) (dotted), and the theoretical boundary term A(r) 
is plotted as a solid line. 



2.3. Intrinsic 'Error' in the PME 

From the analysis of the boundary term in Eqs. (6,15) 
and the numerical calculations performed, we conclude 
that the PME can overestimates the mass by « 20% of 
its true value if the sampling is made only to « li?e of 
the gravitational system. This maximum error is for the 
King's model, while lower values were obtained for more 
concentrated models like Hernquist's or de Vaucouleurs' 
in agreement with the theoretical expectations; see Fig. 1. 

It is important to recall at this point the assumptions 
made in the derivation of the PME: spherical symmetry, 
isotropic velocity distribution and no mass spectrum; im- 
plicitly it also assumes that the particles we observe trace 
matter. If dark matter dominates the outer parts of galaxy 
clusters we expect that sampling the luminous part would 
yield an overestimate in the mass, not an underestimate, 
which can be of « 20%. This particular value depends, 
obviously, on the total mass profile, but can be considered 
typical for profiles suited for galaxy clusters. 

Before we leave this section we address the results 
of Thomas & Couchman (1992) on the PME. Examin- 
ing their formulae B2 we note that these equations cor- 
respond to the case of test particles orbiting a massive 
particle, as they acknowledge, but for non-test particles 
they would require an extra factor of two. It seems, from 
the numbers they give in their Table 5 that this extra fac- 



Hector Aceves & Jaime Perea: On two mass estimators for clusters of galaxies 



tor of two and perhaps a higher constant anisotropy, e.g. 
j3 = 0.75, would help solving the discrepancy they men- 
tion, although it is not clear that one can characterize by 
a constant anisotropy their A^-body clusters. 

3. The Virial Mass Estimator 

In this section we test the virial theorem as a mass es- 
timator. We consider its discrete and continuous form in 
order to understand and clarify some matters regarding 
its use that have recently appeared in the literature (e.g. 
Carlberg et al. 1996). We first outhne the deduction of 
the scalar virial theorem for the sake of completeness and 
future reference (see e.g. HTB). 
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3.1. Discrete Virial Theorem 

Consider a system of N equal-mass particles interacting 
gravitationally. First, we differentiate twice the moment 
of inertia of the system / = (m/2) ^ r^ • r^ with respect 
to time: 



Fig. 3. Differences in virial mass at different radii, both pro- 
jected and non-projected, with respect to the true mass of the 
Dehnen's model are shown. Also, the boundary term A(r) from 
the PME is shown. 
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The second term on the right-hand side is twice the total 
kinetic energy of the system, and the first term is the 
potential energy 
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where r^ == |rj -rj|. 

To derive the common form of the virial theorem we 
require that the time average of the second derivative of 
the moment of inertia vanishes (e.g. Goldstein 1980). In 
this case we have, after averaging over all possible ori- 
entations and assuming an isotropic velocity distribution 
(Limber & Mathews 1960), that 



Mv = 



2G E.<,1/^ 



(17) 



where Rij is the projected separation among pairs of par- 
ticles. This is what we explicitly mean by the virial mass 
estimator, VME. An important observational advantage 
of the VME in determining the total mass of a system is 
that it has virtually no dependence on velocity anisotropy 
for near-spherical systems (e.g. The & White 1986). 

One may also derive a virial estimator by using the dis- 
tances of test particles to the centre of the system (Bahcall 
& Tremaine 1981). Nonetheless, such form of the virial 
theorem is of no interest here since we do not consider 
galaxies in clusters to be properly represented by test par- 
ticles. 



3.1.1. Numerical Test of the VME 

We applied the VME to the numerical models described 
in §2. We proceeded as in the testing of the PME using 
different aperture radii, where we applied the VME to 
the particles inside the corresponding aperture. The mass 
derived using the VME, My, is shown also in Fig. 2 for the 
Dehnen's model. Note that the VME follows more closely 
the projected mass M{R) than the PME over most radii, 
especially for R> R^ . 

Nonetheless, significant deviations occur at _R < _Rc . In 
Fig. 3 we show for the Dehnen's model the difference be- 
tween the VME estimate and the physical mass at differ- 
ent radii [short- dashed), (5v, and that of the non-projected 
virial estimator with respect to M{r) (solid), S^y. The 
maximum differences are 6y ~ 15% and d^y ~ 12%. For 
the King profile we obtained a somewhat higher S^y ~ 
20%, this since it has a shallower profile than the Dehnen's 
one. For the Hernquist's model we obtained lower values 
in general: Sy « 12% and S^y ~ 10%. Also shown in Fig. 
3 is the boundary term (15) that appeared in the PME 
{dot-da.shed). 

The application of the VME to our simulated iV-body 
system does not yield any underestimate of the true mass. 
On the contrary, overestimates for the mass at the inner 
regions of the numerical systems were found. These results 
are contrary to those found by TC in their application of 
the VME, but in the same direction as noted by CYE. 
We will now turn to explain the reason why we have such 
discrepancies in the result from our recent application of 
the VME and M{R) at different radii. 

Here we recall that the basic idea behind the use of 
the virial theorem as a mass estimator is to relate the 
gravitational potential energy and the kinetic energy of 
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the system. If the complete system is in virial equilibrium 
any subsystem has to be also in equilibrium. Computing of 
the kinetic energy poses no problem, but the gravitational 
energy, i.e. the term involving J2 ^/Rij^ has not been cor- 
rectly considered earlier when we computed M\{R). 

Unlike force, the potential energy at a particular point 
f inside a sphere depends on the particles outside this 
radius. In other words, the virial theorem holds at all radii 
if the total potential energy of the particles inside this f 
is accounted for. As said before, the kinetic part of the 
virial theorem does not suffers from this inconvenient. All 
this is important to take into account when mass profiles 
of clusters are to be computed and compared e.g. to mass 
profiles derived from X-ray observations or gravitational 
lensing. 

Therefore, applying the VME to clusters, for which we 
probably do not have a fair value of its true boundary, can 
yield an overestimate of the mass because one ignores the 
contribution to the potential energy of the outer parts of 
the 'total' cluster which may be dark. The neglect of this 
extra potential energy is transformed apparently as a mass 
excess inside the sampled cluster, since the system requires 
this extra mass to be in equilibrium in order to equate 
the corresponding kinetic energy. This is also equivalent 
to having particles of different masses as we increase the 
aperture in the system (heavier particles inside and lighter 
ones outside) . We need to consider a surface pressure term 
in the VME (17) if e.g. the system were confined by an 
external force, but this is not the case here. Note that 
we do not take into account here the probable effect of 
external tidal fields on the mass determination of clusters 
of galaxies. 

In Fig. 4 we show the results of a consistent calcu- 
lation of the potential energy in the VME (eq. 17) that 
supports our previous qualitative explanation. We con- 
sider here only the projected form of this virial estimator 
owed to its use in observational astronomy; an almost per- 
fect agreement was found when we used its non-projected 
form. In dotted line we show the mass estimated using 
the virial theorem by considering particles only within a 
certain 'aperture' radius R; i.e. the summation of inter- 
particle separations in Eq. (17), for each i-th particle, con- 
siders only the j-th particles that their physical position 
is Rj < R. The projected mass of the Dehnen's model is 
shown as a solid line, and the mass estimated using what 
the theory of the virial theorem indicates, Mq, is shown in 
open squares; i.e. the potential energy of the i-th particle 
inside this R is computed by adding the contributions of 
all other particles to ^ l/Rij, including those that have 
Rj > R. 

As observed from Fig. 4, the two above procedures 
yield different results, except for small discrepancies which 
we attribute to numerical artifacts in the construction of 
the iV-body system. The error bars in Fig. 4 correspond to 
Poisson fluctuations oc I/a/TV, and the maximum discrep- 
ancy is « 0.04 which consistent with the numerical noise. 
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Fig. 4. Projected mass profile derived from the virial theo- 
rem (eq. 17) applied to an A'^-body equilibrium and isotropic 
Dehnen system. Ma assumes that the summation of interpar- 
ticle distances is constrained to particles inside a particular 
radius R. The solid line indicates the corresponding projected 
mass, and Mq considers correctly that the particles outside R 
contribute also to the potential energy of particles inside R. 



Similar differences were observed for the other models con- 
sidered here. For systems with less particles the behaviour 
is the same, but of course the errors are larger. 

The physical reason for the overestimation of the mass 
in the application of the PME is also now clear: it does not 
account for the long-range nature of gravity. The product 
v^.Ri in the PME only considers the distance of a galaxy 
to the centre of the mass distribution and does not account 
for its interaction with other particles, as the summation 
involving l/Rij in the VME does. 

The behaviour shown in Fig. 4 supports the idea that 
the VME is an excellent mass estimator. In real appli- 
cations, we can expect the error due only to the VME 
method to be within « 20% for realistic mass profiles 
if sampling is made around the total effective radius of 
the cluster. On other hand, taking Afy (-^^o in Fig- 4) 
as a good estimate of the mass profile its difference with 
A/p may indicate us how far into the equilibrium part of 
a gravitational system are we sampling, and thus lower 
the previous error bound. It might also tell us something 
about its velocity anisotropy. 

In astronomical applications mass estimates are some- 
what unreliable due to the fact that clusters may not be 
in virial equilibrium and/or due to the presence of inter- 
lopers. Also, the determination of the physical extent of 
a cluster is problematic. Theoretically speaking, it is im- 
portant to recall that the VME is not very sensitive to 
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anisotropics (e.g. Aceves & Perea 1998) so its output can 
be more indicative of the equilibrium state of the system 
than of its true mass. On other hand, the use of gravi- 
tational lensing methods, which arc not hindered by the 
requirement of the cluster to be in equilibrium, can in 
conjunction with VME results assess the reliability of the 
latter. However mass determination by lensing methods 
are particularly affected if the cluster under investigation 
is elongated along the linc-of-sight. 

A comment on the TC results for the VME is also 
pertinent here. Their projected mass estimator V2 is again 
only suited when particles (e.g. galaxies) can be treated 
as test particles. If the latter is not true, the application 
of V2 underestimates the actual mass. In fact, when such 
estimator was applied e.g. to our Hernquist A'^-body model 
we obtained V2 « 0.2, a factor 5 lower than the true mass 
which is in agreement with the behaviour observed by TC. 
We now turn to consider the continuous form of the virial 
theorem. 
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3.2. Continuous Virial Theorem 

To derive the scalar virial theorem in its continuous form 
we multiply the Jeans hydrostatic equilibrium equation by 
r , and then integrate over the volume of the gravitational 
system (Binney & Tremaine 1987). Using equation (2) for 
an isotropic velocity system, we have: 



r — — ^ Anr dr 
dr 



dip 2 
rp—— Aiir dr 
dr 



(18) 



Integrating by parts the left hand side and expressing the 
right hand side in terms of the mass interior to r, we have: 



(19) 



r*p(jf. -3 r^paf. dr = - / GrpM{r) dr 



The surface pressure term is iPV = 47rr'^pcr^ (e.g. CYE). 
For future reference, we may estimate the position of the 
maximum of this term for the Dehnen's model, 3PV — 
G'M2ar3(a+6r)(a+r)-VlO : r^ax = (7+\/65)a/8 « 0.65, 
yielding a maximum value of SPVmax ~ 0.04 . 

It is readily verified that the Jeans equation holds at 
every radius, giving e.g. for the Dehnen's model at each 



side of equation (9): —3GM^ar{r 



747r. Thus, ne- 



glecting the 3PV term, which comes from an integration 
by parts, would obviously produce an error in the mass 
estimate when using Eq. (19) (e.g. The & White 1986). 
Now, since the value of the mass M{r) is inside an inte- 
gral we cannot simply subtract the 3PV term point by 
point, or its maximum value, to an estimate of the mass 
using the VME (17) in order to obtain the correct value; 
but see below when this term is expressed in mass units 
(e.g. Girardi et al. 1998). This is because the 3PV term 
is local while the integrals in Eq. (19) and the mass are 
cumulative quantities. 

To quantify the effect of the 3PV term on the mass de- 
termination we should formally solve the integral equation 



Fig. 5. Effect of neglecting the surface pressure {3PV) term 
in the continuous form of the virial theorem. Mnsp indicates 
that no 3PV term is considered in the solution of the integral 
equation (19). Msp takes the referred term into consideration 
when actually solving the integral equation (19) for M{r), and 
M{r) is the theoretical Dehnen mass profile. Neglecting the 
3PV term yields an overestimate of mass at r < i?o while an 
underestimate for r > i?o . In the lower panel the difference 
A/nsp — M{r) is shown for a particular range of values. 



(19). Fortunately, if one does not considers this surface 
pressure term the situation is simple. We have: 



/ { 3r^pa^ - GrpM{r) } dr 



-0, 



(20) 



which we may solve readily, leading to an expression for 
the mass with no surface pressure term as follows: 



MnsAr) = -^ral{r). 



(21) 



From Eq. (21) we can see that the neglect of the 'iPV 
term in (19) does not always overestimates the mass, but 
can even underestimate it in realistic cases due to the ra- 
dial dependence of the velocity dispersion. Indeed, for the 
extreme case of an isothermal sphere, where af. is a con- 
stant, we have an overestimation since Mnsp('') oc r . But, 



e.g. for the Dehnen's model where af. 



at large radii 



the mass will cease to increase settling to a constant value 
at some finite radius. We may estimate the mass interior 
to r = 100, using the expression for a'^ in §2.1.1, giving 
-Mnsp(IOO) « 0.6; i.e. an underestimation of the total 
mass by « 40%. Therefore, in general, an overestimation 
occurs only at the region of the system that may be consid- 
ered nearly isothermal, e.g. r < Re (see Fig. 5), while an 



Hector Aceves & Jaime Perea: On two mass estimators for clusters of galaxies 



underestimation occurs outside of it. The maximum over- 
estimate is approximately given by the numerical value of 

Girardi et al. (1998) have introduced a correction term 
in Eq. (17) when partial sampling of a system is done 
based on an suggestion by CYE; namely, that an overesti- 
mate in mass occurs when applying the VME (eq. 17) due 
to the neglect of the 3PV term in equation (19). Girardi 
et al. have proposed the following formula to correct for 
virial mass estimates when partial sampling is done up to 
a boundary radius b 



Afcv - M\ 



1- 



ATTb^p{b)a^{b) 



J inr'^pdr x a'^(< b) 



(22) 



where cr(< b) refers to the integrated three-dimensional 
velocity dispersion, which we calculate in our numerical 
models as 



f72(< R) = 



N{<R) 



E^: 



2 



where A^(< R) refers to the total number of particles in- 
side radius R and the summation considers only the par- 
ticles inside this R. The correction term in (22) expresses 
the 3PV term in 'mass units' by dividing it by a term 
related to the total kinetic energy of the system up to the 
radius b. 

In Fig. 6 [upper panel) we show the computed pro- 
jected mass profile using the VME by apertures including 
the 3PV correction term (dashed line, Mqv) as in equa- 
tion (22) for our A^-body equilibrium system. The values 
of p{r), (Tr{r), and M{r) were taken here directly from 
the theoretical model, hence representing an ideal obser- 
vational situation, while the integrated velocity dispersion 
was calculated as indicated above. For comparison, the 
mass using the VME when the outer particles are con- 
sidered, Mo, is also shown. In the lower panel the corre- 
sponding differences are displayed. 

From Fig. 6 it follows that both the masses obtained 
from the consistent application of the VME (eq. 17), and 
the correction due to the 3PV term introduced to aper- 
ture values of the VME (eq. 22) provide good estimates 
of the mass profile M{R). However the application of the 
VME to the TV-body system yields a better result, within 
numerical errors, than Eq. (22) . We emphasize that in as- 
tronomical practice the VME does not require to assume 
any particular form of the radial velocity dispersion of 
galaxies, which is needed when constructing the 3PV cor- 
rection in (22). This being perhaps the best feature of the 
virial mass estimator. 

Moreover, in an astronomical situation, in order to 
obtain the 3PV term one has to determine the number 
surface density of galaxies and then apply a deprojection 
integral to obtain the physical number density n(r), or 
use the corresponding p{r) of the fitted model. Here, im- 
plicitly, some confidence is required that the equilibrium 
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Fig. 6. Projected mass profile by the application of the 
VME by apertures, corrected by a surface pressure term for a 
Dehnen's model, Mcv. M{R) denotes the true projected mass 
for the model. In the lower panel, the differences between these 
quantities are shown. The difference M{R) — Mcv cannot be 
accounted by Poisson error bars. 



part of a cluster has been sufficiently sampled to obtain a 
reliable total virial mass M and integrated velocity pro- 
file; as recognized by Girardi et al. (1998). Hence, in this 
respect, the estimation of a mass profile using the the cor- 
rection term in (22) and that on the VME rest on the 
same working conditions, but the latter does not require 
any assumption on ar(r) or its projected version. The 
PME correction term (6) is hindered by the same situa- 
tion. Thus, we are lead to conclude that the VME is a 
reliable estimator of mass and mass profiles, having aside 
the virtue of being straightforward to apply. 

As the PME, the VME can overestimate the luminous 
matter of a cluster if its extent is small in comparison to 
the extent of the mass. In this case one is led to conclude 
that the mass-to-light ratios of the visible parts of clus- 
ters may be importantly overestimated, say by ~ 20%, if 
their outer parts are mainly dark. To conclude this section 
we recall that in deriving equation (17) the assumptions 
were the same as in the derivation of the PME, so these 
equations are valid in this context, as we have verified by 
numerical experiments that have correctly considered the 
potential energy term in Eq. 17. In a future work we will 
consider the estimation of mass profiles of several clusters 
and compare the results with those obtained by other au- 
thors, both in the optical and X-ray region of the spectra. 
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4. Conclusions 

We have investigated the behaviour of the PME and VME 
when applied to a isotropic equilibrium system, and when 
partial sampling of this has been performed. 

We summarize our main conclusions as follows: 

1. The projected mass estimator (PME) and the virial 
mass estimator (VME) yield accurate results for the 
total mass of a gravitational system, provided that the 
whole system is sampled, is in equilibrium and has an 
isotropic velocity dispersion. Moreover, the VME can 
provide an accurate mass profile under the previous 
conditions without requiring the isotropy condition. 

2. The PME overestimates the mass if the sampled region 
is a small portion of the total system. The maximum 
error occurs around the total system's effective radius 
Re, and depends on the density profile. For realistic 
profiles this overestimate can be up to ss 20%. The 
VME yield similar errors when partial sampling of a 
complete system is made. This is because the summa- 
tion ^ ^/Rij in the VME (eq. 17) does not includes all 
members of the system, but only those we can observe. 

3. A recently introduced correction term based on the 
surface pressure term of the continuous virial theorem 
is found to perform, for practical matters, as well as the 
VME under the same kind of hypotheses on the stellar 
system. In applications, e.g. to compute the mass dis- 
tribution of clusters using galaxy data, however, the 
VME should be preferred here since it does not re- 
quire any assumptions on the radial velocity disper- 
sion of galaxies. Discrepancies found in the use of the 
VME, or the PME, in some A^-body simulations arc at- 
tributed to treating particles as test-particles and/or 
to anisotropics present in the numerical clusters. 
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